Entanglement between photons and atoms coupled out from a 
Bose-Einstein-Condensate 



O 

o 

(N 

C 
3 



> 
in 
o 
o 
\o 
o 



X 



T. Gasenzer, D.C. Roberts, and K. Burnett 
Clarendon Laboratory, Department of Physics, University of Oxford, Oxford 0X1 3PU, United Kingdom 

(February 1, 2008) 

We study the limitations to the relative number squeezing between photons and atoms coupled out 
from a homogeneous Bose-Einstein-Condensate. We consider the coupling between the translational 
atomic states by two photon Bragg processes, with one of the photon modes involved in the Bragg 
process in a coherent state, and the other initially unpopulated. We start with an interacting Bose- 
condensate at zero temperature and compute the time evolution for the system. We study the 
squeezing, i.e. the variance of the occupation number difference between the second photon and the 
atomic cm. mode. We discuss how collisions between the atoms and photon rescattering affect the 
degree of squeezing which may be reached in such experiments. 
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Producing entanglement with macroscopic samples 
of atoms has become an increasingly important goal 
for experimental studies with Bose-Einstein-Condensates 
(BECs). Its achievement would, in particular, be of im- 
mense value in quantum information processing (for re- 
views cf. for Heisenberg limited clocks |^| and for 
tests of quantum state reduction 0. Techniques for 
the creation of entangled photon pairs Gj are already 
available for quantum cryptography and teleportation 
P], analogous schemes involving macroscopic samples of 
atoms are still being sought. There are good reasons for 
doing this: whereas photons may be transmitted easily 
over large distances, atoms have the advantage of being 
more easily stored. 

In H we have investigated the possibilities for exciting 
squeezed atoms from condensates using light scattering. 
In this letter we study the entanglement of atoms and 
photons by scattering photons off atoms in a BEC. En- 
tanglement between atoms and photons has been consid- 
ered by others and possible schemes have been proposed, 
In extending the existing work we have studied 
the factors which may limit the amount of entanglement 
achievable, in particular due to interactions between the 
atoms and between atoms and photons. We consider 
weakly interacting atoms with just two internal states in 
a homogeneous BE-condensed state at zero temperature. 
We consider in detail 2-photon Bragg transitions between 
different centre of mass motional states of the atoms in 
the internal ground state, with momenta p, p' and en- 
ergies hujp, TiCjpi |fTof . The intermediate state is taken 
to be an internal excited state. The initial state will be 
assumed to be a product of the weakly interacting con- 
densate ground state and a coherently populated photon 
mode with momentum k^. We then calculate the time 
evolution of the system. Under the assumption that the 
light induced coupling between the atomic modes is weak 
enough to neglect the depletion of the condensate due to 
this interaction. We envisage an experiment, where a cw 
laser with frequency w\ and photon momentum k\ is inci- 



dent on a BE-condcnsate, where photons with frequency 
u>2 are scattered into modes ki . In general, such an exper- 
iment will produce photons scattered in many different 
directions, together with atoms created in modes with the 
corresponding recoil momenta Afc = fci — k%. However, 
we restrict our calculations to a single scattered photon 
mode, assuming that it may be singled out experimen- 
tally. In this way we define the atom-photon mode-pair, 
whose relative number squeezing we consider. By adjust- 
ing the angle between k\ and some fixed Jc2, the frequency 
difference Aw = uj\ — lu2 can be tuned to the difference 
w p i — uj p of the atomic quasiparticle frequencies. This 
allows access of both the quasiparticle and the particle 
regimes of the centre of mass spectrum . 

Before we introduce the model in more detail we give 
a review of relative number squeezing and introduce a 
measure for it. We would like to study the correlations 
between the atoms and the quantized photon mode. The 
two mode squeezing operator has the form 



5(C) =exp +C*ba 



(1) 



Here a and b are the bosonic creation operators in modes 
A and B respectively, and £ is a c-numbcr. The state re- 
sulting from this operator acting on the particle number 
vacuum state would be ideally squeezed, since particles 
are only generated in pairs. It may be written in the form 

oo 

5(C)|0 o ;0 6 ) = sechr^[-e^tanhrf|n o ;n fe ), (2) 

ra=0 

where £ = rexp(i</>), with r and <fi real and positive (cf. 
]l3||). This state is entangled and the distribution func- 
tion over the different occupation numbers n is thermal. 

The squeezing operator (Q) is, up to a diagonal phase 
matrix, equivalent to the time evolution operator describ- 
ing the Bragg transitions we consider, if we restrict our- 
selves to the simplified case where one of the atomic cen- 
tre of mass modes involved in the 2-photon Bragg transi- 
tion is the condensate mode, p = 0, and the other mode is 



1 



p = Ak. Atoms in mode p — Ak and photons in mode k 2 
are thus produced in pairs, and on resonance, Auj = WAfc, 
the time evolution operator is essentially given by ([!]). 

The squeezing is measured in terms of the variance 
of the occupation number difference. For state (0) this 
variance is zero, 

[A(n a - n b )f = ((h a - h b ) 2 ) - ( n a - h b ) 2 

= (An a ) 2 + {An b ) 2 -2((h a h b )) = 0. (3) 



where h a = a) a, h b = tfb, ((h a h b )) = (n a h b ) - 
(h a )(h b ). We define a quantity which describes the rel- 
ative number squeezing between modes a and b: 



f 3 := [A(n a - n b )f / (h a + h b ). 



(4) 



This parameter is equivalent to the one defined in |1J] 
measuring the squeezing considered in |]l5| . For inde- 
pendent classical, i.e. coherently populated states, the 
correlation term (( n a h b )) in ([}]) is zero and the squeezing 
parameter is equal to 1. The thermal state (||) which has 
£3 = is maximally squeezed. Note that this state may 
also be considered as maximally entangled in the sense 
that the relative entropy S re i = Tr a [p a In p a ] is at its 
maximum, where p a = Trb[p a b] is the reduced density 
matrix obtained from the density matrix p ab correspond- 
ing to the state (||). In the context of spin squeezing 
jl4j the parameter (Q) measures the variance of the 3- 
component J3 = ^(h a — h b ) compared to half the total 
value J/2 = \(h a + h b ) of an abstract angular momen- 
tum. It will be instructive to study also the squeezing of 
J 1 = ^(a'b + b^a) and J 2 = i(a^6 — w a) for which we 
define & := [A J,] 2 /( J/2), i = 1,2. 

There are two reasons why perfect squeezing £3 = is 
spoilt in an experiment: Collisional atomic interactions 
and photon (rescattering) processes different from those 
we above restricted our view upon. 

To a first approximation, the ground state of an ho- 
mogeneous interacting BE-condensed atomic gas, i.e. the 
quasiparticle vacuum, may also be written as a state 
of the type given in (|^), with relative particle number 
squeezing between modes of opposite momenta p and —p. 
Thus, even without photon scattering, the squeezing be- 
tween the (empty) mode ki and atomic mode p' = Ak 
is spoilt, with a variance of the relative occupation being 
equal to the variance of the population of mode p' = Ak, 
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'Afc 



)] 2 = [An» Ak ] 2 + 



Moreover, one has to take into account that there are 
other processes of photon induced scattering into and out 
of the condensate mode. Specifically, an atom in mode 
p = —Ak may be transferred into the condensate to- 
gether with a photon scattered from k\ into k 2 , and vice- 
versa. These processes appreciably affect the squeezing 
between p' = Ak and k 2 at longer times as our calcula- 
tions presented below show. 

To study the time evolution of the squeezing we use 
different methods: We may either treat the effective two- 



photon coupling between the atomic cm. states as a per- 
turbation to the Hamiltonian describing the free evolu- 
tion of the interacting condensate and the photons and 
apply perturbative methods. Alternatively we may com- 
pute the time evolution in an exact way within a re- 
stricted Fock space. Due to the Bose-enhancement of 
scattering processes involving the condensate mode, we 
may restrict our view on the atom modes ±Afc. On res- 
onance (Auj = u~JAk) the Hamiltonian restricted to the 
subspacc of these modes and the photon mode &2 is pe- 
riodic in time and we may use the Floquet-approach |l6| ] 
to solve the Schrodinger equation. 

We consider the homogeneous case, where the atomic 
motion is not confined by a trapping potential, but we 
assume that the system is confined in a box of volume 
V, and apply periodic boundary conditions. The homo- 
geneous picture should be a good approximation for the 
central region of the large alkali condensates achieved in 
present day experiments. An extension of our methods 
to the case of trapped condensates is straightforward. 

The effective Hamiltonian describing 2-photon Bragg 
transitions within the different centre of mass motional 
levels of the atom being in the internal ground state 
is obtained by adiabatic elimination of the intermedi- 
ate state. It is written in terms of the operators a and 
b, which satisfy standard bosonic commutation relations 
[Sfe, a\,] — 5kk'> etc., as H — Hq + H af , where 



-Afce 



Jlt + h.c. 



(5) 
(6) 



Here C pq — (h 2 p 2 /2m + n a UQ)6 pq and M pq = n U S p ^ q ) 
are the coefficients of the quadratic Hamiltonian describ- 
ing the motion of the colliding atoms of mass m. The col- 
lisional interaction has been approximated by the contact 
"potential" V(r) = U 5(r), the coupling U being related 
to the atomic scattering length a by Uq = 4Trah 2 /m. 
Any dependence on the polarization of the photons is 
ignored. We have neglected the AC-Stark-shifts, which 
would modify the diagonal terms in the basis of the par- 
ticle operators and could easily be included when nec- 
essary in modelling an experiment. The 2-photon Rabi 
frequency f2 is given in terms of the dipole momentum 
d ge between the internal ground and excited states, of 
the positive frequency part of the electric field describing 
the laser mode I, and the coupling g ge {k2) of the photon 
mode 2 as 



Q = -2\d ge -e l+ \k 1 )\g ge (k 2 )/A. 



(7) 



A = LOkt — ^0 is the detuning of the laser from the internal 
atomic frequency difference luq- 



2 



Hq, which includes interactions in the Bogoliubov ap- 
proximation, may be diagonalized and written in terms 
of quasiparticle operators (3 P given by 



P P = u p bp + Vpb_ p . 



(8) 



The positive coefficients u p , v p of the transformation are 



1- 



-1/2 2 = 
P 



1, with otp = l + y p -y p {2 + 



ip — V .L — u p ; ■ , V p — - U 

Up) 1/2 , Up = \p\/Po, Po = V^ano, n Q being the con- 



densate density. For simplicity we neglect the change in 
Up, v p by the depletion of the condensate due to photon 
scattering. 

The Hamiltonian in terms of quasiparticle operators is 



H = ^ hQpflp'flp + fiuj2a\ 



k 2 a k 2 



P 



Ake 



h.c 



(9) 



with bp — Upfip — Vp/3_ p and quasiparticle frequencies 



tip = {hpl/2m) y p ^2 + y 2 . 



(10) 



Since the recoiling atoms originate predominately from 
the condensate and have a momentum p — Ak, we will 
focus on the squeezing between the photon mode ki and 
the atomic mode Afc. As inital state we choose 1 1 = 
0) = | 0o5 >, defined by: $ p \(/) ;0) = a k2 \ (f> a ;0) = 0, 
6o | 00 ; ) = \/Ao|0o;O), where No is the number of 
atoms in the condensate. 

In the following we present numerical results for the 
squeezing parameters i = 1, 2, 3, for a realistic exam- 
ple of a sodium condensate. We consider specifically a 
homogeneous condensate of N = 10 6 23 Na atoms within 
a volume V = 10 _8 cm 3 , which we choose such that the 
condensate density n — 10 14 cm -3 is of the order which 
is typically observed in experiments. With a = 2.8 nm we 
thus have p = 2.65 • 10 6 m -1 , hp 1 /2m = (2tt) 1.5 kHz. 

We consider the case where the photon recoil energy 
Auj = u>x — ll>2 is tuned to the resonance Atu ~ uj Ak- 
in Fig. [l] we plot the logarithm of the occupation num- 
bers (n% 2 ), (n b Ak ), (n b _ Ak ), for y A k = 5, i.e. for cou- 
pling into the particle regime of the spectrum (|To|), where 



ing the Floquet approach, by restricting the calculation 
to the subspace of these three modes. This is possible 
since on resonance, the interaction Hamiltonian is peri- 
odic in time, with period T = Tr/Co^k- We have chosen 
a two photon Rabi-frequency of = 1 s . Initially the 
atomic occupation is (n b Ak (t = 0) ) = v Ak due to colli- 
sional interactions. Only for t > 0.5 ms the photon and 
atom number occupations are of the same magnitude, 
We have checked that for times shorter 



l k 2 



than ~ 1 ms the curves are equivalent to those obtained 



from a perturbative calculation to order il 2 within the 
whole Fock space of atomic modes. For t > 1 ms a non- 
perturbative rise becomes significant. For (ny_ Afc ) ap- 
proaching iVo = 10 6 however, the calculation reaches its 
limit of physical applicability since we have neglected the 
depletion of Nq. The range of times where the calcula- 
tion is valid may of course be increased by choosing a 
smaller il. 

Fig. ^a shows the squeezing parameter £3. For times 
t < 10 /is it is close to 1, i.e. there is essentially no squeez- 
ing compared to the classical case of independent modes. 
This is due to the fact that the relative atom-photon 
number variance is still dominated by the initial thermal 
variance of the atom mode population, [An b Ak (t = 0)] 2 = 



)({n 



b 

Ak , 



1) 



y Ak( v Ak + 1)' which, in the parti- 



cle regime where v Ak <C 1, resembles the variance of a 
coherent state. (For recoil momenta in the quasiparticle 
regime, Afc < p$ the initial value of £3 would be consid- 
erably larger than 1.) For t > 10 fis the squeezing starts 
to improve due to the atom-photon pair production, but 
only until t ~ 1 ms. The perturbative calculation to 2nd 
order in Q gives a settling of £3 to ~ v Ak ~ 10 -4 (dotted 
line). The exact Floquet calculation within the restricted 
Fock space results in £3 rising again for t > 3 ms (solid 
line). It is interesting to observe what causes the squeez- 
ing to be limited at large times. Neglecting the 2-photon 
Bragg transitions between p = — Ak and p 1 = results in 
good squeezing where the time evolution is described by 
(0) at long times (dashed line). Including the processes 
involving the mode p = — Ak but neglecting the energies 
Aui and &±Ak means that the time evolution is described 
by the operator 

S(C, 1?) = cxp (-C a r b Ak - i3 a r b_ Ak + h.c)) , (11) 

where £ = $ = iflyHVot. In this case the squeezing be- 
tween p' = Ak and k% is reduced already for t > 0.1 ms 
(dash-dotted curve). It is therefore the energy separa- 
tion of the motional levels from the condensate which 
lowers the degree to which the squeezing is suppressed 
and yields the time evolution as shown by the solid line. 

The three mode operator ( p| ) also shows why the 
two-mode squeezing between p = Ak and &2 decreases 
at large t. A change of the occupation number diffcr- 



hp /2m. The curves have been computed us- ence n 



k 2 



Ak 



by minus one is equivalent to the cre- 



ation of a particle in mode p 



-Ak. Thus the vari- 



ance 



'Afc 



)] 2 must be the same as the variance 



(An b _ Ak ) 2 of the population of the p = —Ak mode and 
numerical calculations reflect this. Hence, the creation of 
a perfectly squeezed state w.r.t. modes p' = Ak, ki re- 
quires the mode p = —Afc to be in a number state. This 
is physically realistic only when the mode is empty. How- 
ever the atomic collisions cause any mode with p ^ to 
be occupied to a certain extent, with a thermal variation 
(An b ) 2 — (n b )((n b } + 1), and this limits the squeezing 



3 



right from the beginning. The photon induced transi- 
tions then lead to an increase of ( n b _ Ak ) with time, with 
(An b _ Ak ) 2 remaining a thermal variance. For the colli- 
sionless case (a = 0) a similar squeezing deterioration at 
large times has been found in B. 

We still note that the squeezing of the difference of the 
quasiparticle occupation number of the atomic mode Ak 
and the particle occupation of the photon mode ki is at 
maximum at time zero. But due to the population of 
the quasiparticle mode —Ak, the time evolution of the 
quasiparticle squeezing is similiar to that of the particle 
squeezing for times t > 3 ms. 
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FIG. 1. Time dependence of the mean particle occupation 
numbers of the scattered photon mode, (n k2 ), and of the 
modes of the recoiling atoms, (n b Ak ), (n_ A/ . }. 



of atoms and photons in the modes Ak and k-i being 
rescattered before leaving the condensate as r x = a x nol, 
where a x is the total cross section for atoms, cr a t = 87ra 2 , 
or photons, er p h = (3/27r)(A7/A) 2 , and / the width of 
the sample. With 7 = 0.6 • 10 8 s _1 being the natural 
linewidth of the sodium A = 589 nm line, and a typi- 
cal detuning of A = 1 GHz we obtain for our example 
condensate r at = 42% and r p h = 3%. This shows that 
although the photon rescattering plays a relatively small 
role, there will be further limitations of squeezing arising 
from loss in the scattered atomic mode, which we have 
not included into our calculations so far. Despite this our 
results show that there is a substantial initial limitation 
of squeezing due to the collisional interactions between 
the atoms, which are naturally present in the gas before 
the laser is switched on. 
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FIG. 2. (a) Time dependence of the squeezing parame- 
ter £3: exact calculation (solid line), perturbative calculation 
(dotted line), energy differences neglected (dash-dotted line), 
secondary recoiling mode p = — Ak neglected (dashed line), 
(b) Time dependence of the squeezing parameters £i. 



In the context of spin squeezing, we found analytically 
that within the restricted Fock space of the modes ±Afc, 
&2, we have ( J%) = ( J2) = for all times. Thus the 
mean spin vector ( J ) points in 3-direction. Fig. |p3 then 
shows, that although the squeezing of all three spin com- 
ponents becomes worse for large times compared to the 
case of coherent modes, there is a strong relative squeez- 
ing of several orders of magnitude of the 3-component 
compared to the 1- and 2-componcnts which are of equal 
size. 

It might be argued that recollision processes of pho- 
tons created in mode k-i may be avoided by choosing 
an optically thin BE-condensate, which the photons are 
leaving before they may recollide. Such a scheme has 
been proposed in This can help to avoid the squeez- 
ing being spoilt at large times. We estimate the fraction 
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